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I. INTRODUCTION 

Since the formulation of General Relativity, about one hundred years ago, many alternative geometric theories have 
been proposed in order to explain gravitation phenomena (e.g. [iHa]). Some of these touch one important aspect 
of General Relativity: the conservation law. To our knowledge, the first non-conservative theory of gravity was the 
steady-state model [1, , following some ideas already presented in the end of the forties by Jordan [llj . In the 
beginning of the seventies, Rastall proposed one new version of a non-conservative theory of gravity, following the 
remark that the conservation law T^'';^ = may not hold true in a curved space-time [E, Hence, he argued that 
new gravitational equations can be obtained considering a modification of the conservation law such that 

T^-",^ - kR'', (1) 

where T^'' is the energy-momentum tensor, k is a coupling constant and R is the Ricci scalar curvature. Hence, in 
the weak field limit, the usual expressions are preserved. Since, generally, the Ricci scalar curvature is connected with 
the trace of the energy-momentum tensor, Eq. ([T]) can be re-written as 

T^'',^ = RT''' , (2) 

where k is a new constant and T is the trace of the energy-momentum tensor. It is curious to remark that the 
phenomenon of particle creation in cosmology (l^ - [T^ also leads to a violation of the classical conservation laws and, 
in this sense, Rastall's idea may be viewed as a kind of classical formulation of that quantum phenomenon, since the 
violation of the energy-momentum conservation is connected with the curvature. 

More in detail, Rastall's modification to Einstein equations take the following form (c = 1 units): 

Ri^f - ]^gtivR = SttG - ^^^Y^g^^T^ , (3) 

T^^,. = ^7^" > (4) 

where 7 is a parameter (the choice 7 = 1 restores General Relativity) . Note that it seems possible to have a Lagrangian 
formulation from which the above equations are deduced jT^ . Since for a radiative fluid T = 0, implying i? = 0, we 
can expect that the cosmological evolution during the radiative phase is the same as in the standard cosmological 
scenario. At same time, a single fluid inflationary model, described by a cosmological constant, is the same as it would 
be in the General Relativity case. Hence, Rastall cosmologies may have an irnportant departure from the standard 
cosmological model from the beginning of the matter dominated phase on |la-[l9j . 

II. THE MODEL 

In order to construct a Rastall cosmology which accounts for the matter dominated era and the present phase of 
acceleration of the universe, let us consider a two-fluid model. The first component is a pressure-less matter (i.e. 
cold dark matter plus baryons) with density p™, while the second one obeys the vacuum equation of state px = —px- 
A subscript (or superscript) m shall denote quantities related to the matter component whereas a subscript (or 
superscript) x shall refer to dark energy quantities. 



Equations ^ and Q then become 



7-1 



(5) 



(rr + TD;, = ^(Tm + TxV . (6) 

The latter equation fix the divergence of the total energy-momentum tensor, while extra assumptions are necessary in 
order to fix (T^^'');^ and (T'^"^);;j (similarly to General Relativity, where the Bianchi identities only constrain the total 
energy-momentum tensor). In general, for arbitrary linear combinations of r,„ and Tx such that Eq. ([6]) is preserved, 
we introduce the arbitrary real parameters r]x and 77™ as follows: 

T^^r^ = ^ (^™^™ + n.Tx)''^ , = ^ [(1 - nrn)T„^ + (1 - V.)Txr , (7) 



Note that, for 7 = 1, we recover the standard (General Relativity) case. 
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Together with the flat Robertson- Walker metric, 

ds^ = dt^ — a'^{t)Sijdx^dx^ 
the Einstein and conservation equations (O , ^ take the form 

SttG 



iHprn 
Px 
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7 - 


1 


2 




7 - 


1 



(3 - 2-f)p^ + ^-7— ^Pm 



[(1 - r]m)Pm + 4(1 " ■llx)Px] 



Combining the second equation with the third, one obtains 



p,i 



3Hp„i — 



7-1 



Vm + 2(7 - l){r]x - r]m) 
1 + 2(7-1)^-1) 



Px = 



7 



1 



1 



Thus, we have the following solutions for the densities: 



Pm ~ PmOfl 
P 



2 l + 2(7-l)(77, 

rim + 2(7 - l)(r?a: - Vm) 



1 + 2(7 - 1)(77, -l)-2^[rim + 2(7 - 1)(77, - 77™)] 



Pa: 



7-1 



1 - Vn 



3-27 



2 l + 2(7-l)(77, -1)' 



(8) 

(9) 
(10) 
(11) 

(12) 

(13) 
(14) 



where p is an integration constant and we have introduced uJe as an effective equation of state. Inserting Eq. (jl4p in 
the modified Friedman equation we obtain 



SttG 



1 



(7 - 1) [(3 - 27) (1 - r;™) - 1 - 2(7 - l){r,, - 1)] 



2[l + 2(7-l)(,7,-l)] 



(15) 



which would describe the same background evolution as in the ACDM model, with now p playing the role of an 
effective cosmological constant, if matter had not now a different evolution, depending on 7, r]m and rj^- 

Equation (|13p predicts a deviation, given by the parameter cjg, from the usual dependence of the matter component 
with respect to the scale factor. Indeed, in the ACDM model: ujg = 0. In order to constrain Wg, we use supernova 
typ e la data, Union2 sample [20| . H{z) data [2l|, ll^j, first CMB acoustic peak da ta g iven in terms of the R factor 
[23} . and BAO data [2^, 'M]. We perform a bayesian statistical analysis (as done in (l7j ) for the model defined by the 
cosmological term plus a fluid with — ^ePm- Three free parameters are consider: the matter density, the Hubble 
constant, and Wg. Marginalizing over the matter density and the Hubble constant, we find, at 2cr: uj^ = 0.039j:o:o80- 
The PDF for lu^ is displayed in Fig. [TJ 
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FIG. 1. Probability density function (PDF) of oje- 

Another estimation can be obtained by using structure formation. Following (25j , we see that the modification of the 
equation of state induced by Rastall's theory at the background level, from the point of view of the general relativity 
framework, i.e. Pm = '^ePm, is maintained for linear perturbations, i.e. 5pm ~ uJeSpm (we use the approximate sign 
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since in [25| the model considered is a single-fluid one). The cold dark matter scenario implies We = 0, while the warm 
dark matter scenario could lead to (we admit that the hot dark matter scenario is excluded by observation). 
Suggested by these remarks, we fix = from now on. 

It is remarkable that even if 7 is significantly different from one, the cosmological background of Rastall's theory 
may lead to the same results of the ACDM model. To this end, we set uje — 0, whose single 7-independent solution is 

Vx = r]m = ■ (16) 

The above values for rj^ and rjm lead to a particular but interesting case of Rastall's theory. The above setting 
guarantees that the matter energy-momentum tensor is conserved in the usual sense, i.e. T^':'^ = . 

In fact, with the above hypothesis, Eq. ([T5|) reduces to the corresponding one in the ACDM model, without any 
explicit dependence on 7. 



III. SMALL PERTURBATIONS 



We turn now to the fate of linear perturbations. We will work initially with the synchronous coordinate condition, 
since it leads directly to an expression for the density contrast of matter perturbations (26l.[27t. It is more convenient, 
in this case, to use the field equations under the form. 

The synchronous coordinate condition implies that 

Oiii' = gtiu + hfj_^ , /i^o = , (18) 

where g^^, is the total (perturbed) metric, (7^^ is the background metric and h^^, are fluctuations over the latter. 
The fluids are represented by the following energy-momentum tensors: 

rr = (p. + P.XK - P.g^" , = (p™ + PrnKn< - Pmg^" ■ (19) 

Later, we shall impose pm = and Px = —px- The choice of the synchronous coordinate condition implies Su'^ = 0, 
but (5u* is a dynamical variable. There are in fact two four-velocities, one associated to the x component, the other 
to the matter one. Using this formalism, one obtains the following set of equations for perturbations: 

h + 2-h^ SnGpxh + 3(2 - 7)0;^,] 4 + 87rGp™7<5m , (20) 
a 

(21) 

5px + (1 +ujx)Px[Q - 77 )+3-(l +L0x)5px ^ ^ p. ^ [(1 - 'iu)x)5px +(5/9^1 , (22) 
\ I I a 2 



(1 + LUx)i3x5v'x + (1 + ^x)PxSii'x + 5-(l + uJx)Px5v^x + ^d'^px = --—^ [(1 - 'i^x)d'Spx + 5'<5p™] . (23) 



In these expressions we have used the deflnitions 

h = , ujx = — , dx = , dm = , O = d^dul , (24) 

a Px Px Pm 

where di denotes derivative with respect to the co- moving spatial coordinates. Now, let us impose ujx — —1. We 
obtain, from Eqs. and ([231), that 



6Px = ^ {^5px + 5prn) , = {^d'Spx + d'Spm) , (25) 



and both these equations lead to the relation 



'^^^ = 2(1-27) '^^'" ■ (2^) 
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Using now Eq. pif) one can rewrite Eq. ([20)) as 



a . 



5m + 2-5m - ^TTGpmSm = . (27) 

a 

But this is the same equation for matter perturbation as in the ACDM model! Actually, there is a difference: now 
there are perturbations in the dark energy term px- In fact, using Eqs. ()13p and ()14|) with ?7,„ — rjx = ^ together with 
we find 

2 -^a^ + ^ ■ 

PmO ^ 

In the remote future (i.e. for a — > oo), 5x must become negligible, and a complete equivalence between the Rastall 
and the ACDM model is expected both for the background expansion and the linear perturbations evolution. On the 
other hand, in the remote past (i.e. for a — 0), we have 5x ~ and therefore dark energy may in principle cluster. 
This may have consequences for the structure formation process because we should expect an amount of clustered 
dark energy to be present in virialized systems, like halos of galaxies and clusters of galaxies. However, as we have 
shown, such consequences seem not to appear for linear perturbations. 



A. Newtonian gauge 



A possibility to discriminate between the two models could be a different integrated Sachs- Wolfe effect signal [28| . 
To see if this is the case, we adopt now the conformal-Newtonian one. Consider the following perturbations to the 
Robertson- Walker metric expressed in the conformal time 77 : 



ds^ = a^(?7)(l + 2<^)dif - a^{r^){l - 2<^)Sijdx'dx' , 



(29) 



where $(77, x^) is the gravitational potential. For each fluid component we calculate the perturbations to the continuity 
equation and to the Euler one. Following [2^, for the case of matter we have 



Sp[^ + 3—5p,n + apmdiSu'„, - 3p™$' = , ^ (a'^pmSul^)' + pmd"-^ = , 



(30) 



where the prime denotes derivation with respect to the conformal time. For the x component, we have the correction 
proportional to T^^ and the equations are 



a! 7 — 1 

ip'x + i—{5px + Spx) - 3pxXl + Wx)<^>' + apx{l + Wx)dM'x = [Sp'm + Sp'x - Wx] 

a 2 



1 



4 i^a^Pxil + wx)6ui] + d'Spx + Px{l + Wx)d'<i> 



7-1 



9* [Sprn + Spx - 3Spx 



For the Einstein equations, we have 



A$ - 3— (—$ + $'] = 4:TTGa^ 



$"-H3-$' 



a \ a 
a 



Spm + Spx 



Spx 



7-1 
2 

7-1 



{Sprn + Spx - SSpx) 
{Sp„i + Spx - ^Spx) 



For px = —px, Eqs. (PT|) and ([5^ become 

7-1 



{Sp',^ + 4Sp'x 



d'Spx 



7-1 



d' {Spm + "iSpx) , 



(31) 
(32) 

(33) 
(34) 

(35) 



and therefore we have 



dpx 



7-1 
2(3-27) 



Spm , 



(36) 
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i.e. the perturbation in the x fluid is proportional to the one for matter, the same relation (j26|) we found using the 
synchronous coordinate condition. Combining Eqs. ([55]) and with Eq. ([55]) we obtain 

(37) 
(38) 



A$ - 3— — $ + $' = inGa^Spn 
a \ a 



$ 



These equations are identical to the corresponding ones for the ACDM model (note that the background is also the 
same). In particular, the first one confirms that matter perturbations are not affected by the x fluid (with lo^ = ^l); 
even if the latter also agglomerates. 



IV. NON-LINEAR REGIME 



We follow '3Cf| for the treatment of the second order regime of cosmological scalar perturbations. In particular, 
consider Eqs. (4.15)-(4.17) of [131 for the x component of our model in the Newtonian gauge (i.e. zero metric coefficient 
B), with zero anisotropic pressure and Px — —Px, ^Px = —Spx (we assume this relation at all orders): 

(^)7^o = Sp^^^ , (^)7^, = , (^)7^o = , (^T^^ = ^Sp(pSl , (39) 

where the subscript (■) refers to the perturbation order. Moreover, the second-order trace has the following form: 



(2)T. =(2) T'xo +^'^ - Sp(^^ - S'JP?) = 4<5p?) . 
On the other hand, for the matter component we have p,„ = 0. Hence Eqs. (4.15), (4.17) of [sOl become 

where the v'"^^ represents the velocity perturbation at first order. Hence, for the second-order trace we get 

(2)^ _ (2)^0 , (2)^fe _ e (2) 

Now consider the equation 

7-1 



at second order: 



f) (2)71^ I (O)pA' (2)tiP I (l)pA' (l)7^P I (2)pA' {0}rpp 

M X ]j ' /ip X u ' pp X u ' pp X I' 

_ (0)pp (2)^p _ (l)pp (l)^p _ (2)pp (0)™ 

/.il/ X p f.LU X p pU X p 



7-1 



(2)rp . (2)rp 



Making = 0, and taking into account the non-vanishing values of the perturbations, we obtain 
3^ (2)7^0 , (0)pp (2)7^0 , (i)pp (1)7,0 , (2)pp (0)7,0 

•^0 2; ^ pO 2; ^ pO 2; ^ pO ^ K 

7 - 1 [(2)7, 



_ (0)pP (2)7np _ (l)pp (l)7,p _ (2)pP (0)7np ^ 

-^pO xp p-a xp -^pO xp 2 



i^)Tx 



On the other hand: 



(0)pP (2)^p _ (0)pO (2)^0 , (0)pj (2)rj., 
pO ^x n — no J- r n ~r i ^x 



^x p — ^00 ^ a; ""^ jO 

t' r 



a 
a 



X I 

/ 



and 



(40) 



(41) 



(42) 



(43) 



(44) 



(45) 



(46) 



(47) 
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Hence, the second term in the first line of Eq. (|45|) . cancels with the first term of the second line. In the same way, 
the third term and fourth term of the first line cancel with the second and third of the second line. We end up with 



(2)' 



7-1 



45^(2) +sp. 



(2)' 



(2) 



2(3-27)^™ ' 



(48) 



and we recover, up to an integration constant, Eq. ([55]) . now at the second order. Therefore, it seems that in the 
Rastall cosmological models here investigated, the evolution of perturbations in the x fluid mimic those in the matter 
up to the second-order perturbation regime. 

The left-hand side of Eq. (jS]) is the same as in Einstein's theory, for the background and for the perturbations 
evolution. As for the right hand side, using Eq. (4.16) of [30] for matter {pm = 0), i.e. 



(2)7-0 . 



(49) 



we find 



(2)^0^ _ (7_1) (2)r ^ ^^(2) ^ 2p„4%(l)'= 



(50) 



This is the same relation that one can find in the usual ACDM model. Following the same steps, we have 



(2)yi 



(2)rp0 _ {2)rp0 _ 

I ~ -'-ml " 



(7-1) 



(1)^,(1) 
j 



(51) 
(52) 



All the relations are the same as in the ACDM model. Therefore, even at the second order, considering scalar 
perturbations only, the evolution of matter perturbations is not affected by the agglomeration of the x component. 

However, in the full non-linear regime one should expect some differences between Rastall's model and the ACDM 
one. This because we have now two components agglomerating and they would gravitationally interact affecting in 
some manner, hopefully detectable, the growth of structure. Such analysis is beyond the scope of the present paper, 
nevertheless we can argue some hints by considering the simple case of a spherical top-hat collapse. Let us follow the 
growth of an inhomogeneity in the past, when the universe behaved as an Einstein-de Sitter one since the effects of a 
cosmological constant were safely negligible. Following [3l| , the growth of a spherical homogeneous region is described 
by the equation 



Hf 



np{u)- + i-np{u) 



(53) 



where ai is the scale factor computed at some time ti at which the collapse begins and fip = 1 + 5 is the density 
parameter of the collapsing region, determined by the density contrast of the fluid there contained. Here comes the 
relevant difference. For the ACDM model, the only collapsing component is matter and, therefore, ftp = 1 + 6,n- On 
the other hand, in the Rastall models we presented, an unknown fluid x with a vacuum equation of state also can 
agglomerate. For this reason, the density parameter of the collapsing region assumes now the form 



nr. 



Spx 



(54) 



Considering now Eq. (|14p with rjm = rjx = 0, and neglecting the contribution of the constant p, we can write 



2(3-27) 7-1 . 
SZ„ = 1 -I z z dm + :^c).T 



5-37 



5-37 



(55) 



Equation ([551) is assumed to hold true for both the models and it is easy to solve, giving the size of the perturbation, 
once we know Up{ti). For the ACDM we have that rip(ti) = l + 6„i{ti), whereas in the Rastall- type cosmological model 
we have one degree of freedom more than the ACDM one, because of the initial condition SxiU). So, in principle, one 
can construct the same collapsing history as in the ACDM one but with a different amount of matter, depending on 
the value of 7. 
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V. DISCUSSION AND CONCLUSIONS 



We have shown that there is a subset of cosmological scenarios based on Rastall's energy-momentum tensor non- 
conservation that are oquivalcint to the ACDM cosmology, except for one aspect: vacuum energy may agglomerate. 
Two effects could allow to discriminate between the two models: non-linear effects in the matter power spectrum and 
the transfer function for cosmological perturbations. Vacuum energy is negligible in the past, hence the impact of its 
fluctuation in the evolution of the perturbations in the primordial periods of the evolution of the universe may not 
be relevant. However, according to the scenario described above, dark energy must be present in virialized systems, 
like galaxies and cluster of galaxies, and the effect of agglomeration of dark energy must be relevant at this level. 
Therefore, Rastall's cosmology and the ACDM model seem to be distinguishable only at the non-linear regime of the 
evolution of cosmic perturbations. Though we have given some hints of the latter possibility addressing the simple 
case of a spherical top-hat collapse, a deeper analysis shall be the subject of future investigation. 
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